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Abstract 

It is widely accepted that the prompt transient signal in the 10 keV - 10 GeV 
band from gamma-ray bursts (GRBs) arises from multiple shocks internal to the 
ultra-relativistic expansion. The detailed understanding of the dissipation and ac- 
companying acceleration at these shocks is a currently topical subject. This paper 
explores the relationship between GRB prompt emission spectra and the electron (or 
ion) acceleration properties at the relativistic shocks that pertain to GRB models. 
The focus is on the array of possible high-energy power-law indices in accelerated 
populations, highlighting how spectra above 1 MeV can probe the field obliquity in 
GRB internal shocks, and the character of hydromagnetic turbulence in their envi- 
rons. It is emphasized that diffusive shock acceleration theory generates no canonical 
spectrum at relativistic MHD discontinuities. This diversity is commensurate with 
the significant range of spectral indices discerned in prompt burst emission. Such 
system diagnostics are now being enhanced by the broadband spectral coverage of 
bursts by the Fermi Gamma-Ray Space Telescope; while the Gamma-Ray Burst 
Monitor (GBM) provides key diagnostics on the lower energy portions of the parti- 
cle population, the focus here is on constraints in the non-thermal, power-law regime 
of the particle distribution that are provided by the Large Area Telescope (LAT). 
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1 Introduction 



Despite rapid evolution of the understanding of gamma-ray bursts over the 
last decade since the first redshift determinations established a cosmological 
distance scale, there is still much to be learned about associations, progenitors, 
and the radiation processes active in both the prompt and afterglow emission 
regions. The most popular burst paradigm for the genesis of the prompt burst 
emission is of radiative dissipation at internal shocks that accelerate particles 
(Rees & Meszaros 1992; Piran 1999; Meszaros 2002). Within this scenario, 
it is of great interest to understand what physical conditions in the shocked 
environs can elicit the observed high energy indices and the spectral structure 
around the MeV-band peak. It is clear that the measurement of the high 
energy spectral index provides a key constraint on the interpretation of the 
electron acceleration process. To provides insights into these conditions, one 
turns to models of diffusive shock acceleration in relativistic systems. 

A core property of acceleration at the relativistic shocks that are presumed 
to seed prompt GRB radiation is that the distribution functions /(p) are 
inherently anisotropic. This renders the power-law indices and other distri- 
bution characteristics sensitive to directional influences such as the magnetic 
field orientation, and the nature of MHD turbulence that often propagates 
along the field lines. Consequently, familiar results from relativistic shock ac- 
celeration theory such as the so-called canonical a = 2.23 power-law index 
(e.g. Kirk et al. 2000) are of fairly limited applicability, though they do provide 
useful insights. This diversity is fortunate since the GRB database so far has 
exhibited a substantial range of spectral indices, so any universal signature in 
the acceleration predictions would be unnecessarily limiting. 

This paper explores some of the features of diffusive shock acceleration using 
results from a test particle Monte Carlo simulation, and addresses probes of the 
theoretical parameter space imposed by extant GRB observations above the 
MeV spectral break. The Monte Carlo approach (Ellison, Jones & Reynolds 
1990; EUison & Double 2004; Niemiec & Ostrowski 2004; Stecker, et al. 2007) 
is one of several major techniques devised to model particle acceleration at 
shocks; others include semi-analytic solutions of the diffusion-convection equa- 
tion (Kirk & Heavens 1989; Kirk et al. 2000), and particle- in-ccll (PIC) full 
plasma simulations (Hoshino, et al. 1992; Nishikawa, et al. 2005; Medvedev, et 
al. 2005; Spitkovsky 2008). Each has its merits and hmitations. Tractabihty 
of the analytic approaches generally restricts solution to power-law regimes 
for the /(p) distributions. PIC codes are rich in their information on shock- 
layer electrodynamics and turbulence. To interface with GRB data, a broad 
dynamic range in momenta is desirable, and this is the natural niche of Monte 
Carlo simulation techniques, the focus of this paper. 
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Useful diagnostics (Baring & Braby 2004; Baring 2009) on /(p) have al- 
ready been enabled by data from the BATSE and EGRET instruments on the 
Compton Gamma-Ray Observatory (CGRO) for a few bright bursts. Signifi- 
cant advances are anticipated in the understanding of such constraints in the 
next few years, afforded by the broad spectral coverage and sensitivity of the 
GBM and LAT experiments on NASA's Fermi Gamma-Ray Space Telescope; 
this is already being realized. 



2 Diffusive Acceleration at Relativistic Shocks 

The simulation used here to model diffusive acceleration in relativistic pla- 
nar shocks is a kinematic Monte Carlo technique that has been employed ex- 
tensively in supernova remnant and heliospheric contexts, and is described in 
detail in numerous papers (EUison, Jones & Reynolds 1990, hereafter EJR90; 
Jones & Ellison 1991; Ellison, et al. 1995; Ellison & Double 2004; Summer- 
lin & Baring 2006; Baring 2009). It is conceptually similar to Bell's (1978) 
test particle approach to diffusive shock acceleration. Test particles that are 
injected upstream gyrate in laminar electromagnetic fields, their trajectories 
being governed by a relativistic Lorentz force equation in the frame of the 
shock. In general, the fluid frame magnetic field is inclined at an angle ©Bfi 
to the shock normal. Because the shock is moving with a velocity u(x) relative 
to the plasma rest frame, there is, in general, a u x B electric field in addition 
to the bulk magnetic field. Particle interactions with Alfven wave and other 
hydromagnetic turbulence are modeled by using a phenomenological scatter- 
ing of the charges in the rest frame of the plasma. The scattering precipitates 
spatial diffusion of particles along magnetic field lines, and to a varying extent, 
across them as well. The scatterings are also assumed to be quasi-elastic, an 
idealization that is usually valid because in most astrophysical systems the 
flow speed far exceeds the Alfven speed, and contributions from stochastic 
second-order Fermi acceleration are small. The diffusion permits a minority of 
particles to transit the shock plane numerous times, gaining energy with each 
crossing via the shock drift and first-order Fermi processes. 

A continuum of scattering angles, between large-angle or small-angle cases, 
can be modeled by the simulation. In the local fluid frame, the time, 6tf , be- 
tween scatterings is coupled (EJR90) to the mean free path, A , and the maxi- 
mum scattering (i.e. momentum deflection) angle, 6'scatt via 5t/ ~ -^^scatt/ (6"^) 
for particles of speed v ^ c. Usually A is assumed to be proportional to a 
power of the particle momentum p (see EJR90 and Giacalone, et al. 1992) for 
microphysical justiflcations for this choice), and for simplicity it is presumed to 
scale as the particle gyroradius, , i.e. A = rjrg oc p . The parameter rj in the 
model is a measure of the level of turbulence present in the system, coupling di- 
rectly to the amount of cross-field diffusion, such that rj — 1 corresponds to the 
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isotropic Bohm diffusion limit, where the field fluctuations satisfy SB/B ~ 1 . 
In kinetic theory, 77 couples the parallel ( = Xv/3) and perpendicular ( k,j_ ) 
spatial difi^usion coefficients via the relation k±/k\\ = 1/(1 + rj"^) (Forman, et 
al. 1974; Ellison, et al. 1995). In parallel shocks, where the B field is directed 
along the shock normal ( ©Bfi — 0), V has only limited impact on the resulting 
energy spectrum, principally determining the diffusive spatial scale normal to 
the shock. However, in oblique rclativistic shocks where Obh > , the diffusive 
transport of particles across the field (and hence across the shock) becomes 
critical to retention of them in the acceleration process. Accordingly, for such 
systems, the interplay between the field angle and the value of rj controls 
the spectral index of the particle distribution (Ellison & Double 2004; Baring 
2004), a feature that is central to the interpretation of GRB spectra below. 

The test particle assumption adopted in this paper is appropriate as long 
as the energy density Ua of accelerated particles is much less than that of 
the thermal gas Ug . This is generally the case for the simulation results pre- 
sented here. For the most energetic particles to establish Ucr ^ Ug demands 
very fiat distributions, namely a < 2 if dN/dp = 47rp^/(p) oc p~'^ . It can 
be inferred from the distributions in Fig. 1 that those cases that satisfy this 
condition are inefficient at injecting from thermal energies, so that they gen- 
erate Ug/Ucr <^ 1 . These contrast non- linear acceleration scenarios where the 
non-thermal population modifies the Rankine-Hugoniot MHD structure of the 
shock because ^ Ug , thereby inducing spectral concavity in the energetic 
tail. Such non-linear cases have been explored extensively for the contexts of 
galactic cosmic rays produced in supernova remnant shells (e.g. Ellison, et 
al. 2000) and the Earth's bow shock (Ellison, Mobius & Paschmann 1990, 
hereafter EMP90) Generally, as is apparent in Fig. 2 below, GRB spectra are 
too steep to sample the non-linear acceleration regime. Moreover, in the case 
of GRB 080916c, the power-law character above 1 McV is well established 
(Abdo, et al. 2009a): no evidence of any spectral concavity can be discerned. 
While some models of GRB dissipation in the literature assume that large 
fractions (i.e. > 30 %) of the total lepton energy density reside in the accel- 
erated electron population (perhaps adopting truncated power-laws), there is 
no theoretical mandate for such from the perspective of diffusive shock accel- 
eration theory. This assertion applies not only to results from Monte Carlo 
codes such as are presented here, but also to distributions generated in PIC 
simulations (e.g. Spitkovsky 2008). 

The mildly-relativistic shock regime u < c forms the focus of this exposi- 
tion, since this is germane to internal GRB shock models that invoke shock 

formation via the collision of two ultra-relativistic shells. The discussion will 
first outline the core characteristics of the diffusive acceleration process, before 
moving onto constraints on theoretical parameters imposed by the observed 
spectral indices of prompt emission in bright gamma-ray bursts. 
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2. 1 Acceleration Characteristics for Relativistic Shocks 

As mentioned above, the key property of difTusive acceleration at relativistic 
shocks that distinguishes them from their non-relativistic counterparts is their 
intrinsic anisotropy. This is driven by the powerful convcctive influence that 
enables efficient loss of particles away from and downstream of the shock. 
The result of this loss is a general difficulty in generating flat distributions of 
shock-accelerated particles, particularly for so-called superluminal (oblique) 
relativistic discontinuities. Here, the array of possible distribution indices a 
is highlighted, spawned by the sensitivity of both the energization in, and 
escape from, the shock layer, to (i) the size of the momentum deflection angle 
^scatt 5 (ii) the frequency or relative mean free path A/r^ of scatterings, and 
(iii) the upstream field obliquity ©sfi , a quantity derived from the global 
MHD structure of the shock. In this paper, the focus will be on the latter two 
influences. 

Before investigating them, it is appropriate to mention the flrst effect, which 
was originally identifled for relativistic shocks in EJR90. When the diffusion 
in the shock layer samples large held fluctuations SB/B ~ 1 , it corresponds 
to large momentum deflections, delineating the regime of large angle scatter- 
ing (LAS) with A/Ti < ^scatt ^ tt , where Fi is the upstream flow's incoming 
Lorentz factor. Such large deflections produce huge gains in particle energy, 
of the order of in a single scattering and therefore also in successive shock 
crossings. These gains are kinematic in origin, and are akin to those in inverse 
Compton scattering. The result is an acceleration distribution dN/dp that is 
highly structured and much flatter on average than p~'^ , flrst noted by EJR90. 
The bumpy structure is kinematic in origin and becomes more pronounced for 
large Fi (Baring 2004; Ellison k Double 2004; Stecker, et al. 2007, hereafter 
SBS07; Baring 2009). For ultra-relativistic shocks, when p » mc, the bumps 
asymptotically relax to form a power-law distribution dN/ dp oc , as is 
mandated by the concomitant lack of a momentum scale, with an index in the 
range of (j ~ 1.6 (SBS07). From the plasma physics perspective, magnetic tur- 
bulence in relativistic shocks could easily be sufficient to effect scatterings on 
intermediate to large angular scales ^g^att ^ 1 /^i , ^ proposition that becomes 
more enticing for ultrarelativistic shocks. 

Particle distributions for ^^scatt ^ V^i are much smoother in appearance, 
and often necessarily steeper, at least for superluminal regimes. Intermediate 
scattering angles 9 

scatt 1/^1 gGHGrate smooth distributions (SBS07, Baring 
2009), much like those for (but flatter than) the more familiar small angle 
scattering (SAS, often called pitch angle diffusion, PAD) which will consti- 
tute the regime of focus hereafter. This regime has spawned the often cited 
asymptotic, ultrarelativistic index of cr = 2.23 for dN/dp oc p^'^ (Kirk et al. 
2000; see also Bednarz & Ostrowski 1998; Baring 1999). This special result is 
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realized only for parallel shocks with 0Bfi = 0° in the mathematical limit of 
small angle scattering ^^scatt ^ ^/^i ! where the particle momentum is stochas- 
tically deflected on arbitrarily small angular (and therefore temporal) scales. 
In such cases, particles diffuse in the region upstream of the shock only until 
their velocity's angle to the shock normal exceeds around 1 /Fi , after which 
they are rapidly swept downstream of the shock. The lower kinematic energy 
gains in shock transits more than compensate for the accompanying slightly 
higher shock-layer retention rates, producing a steeper distribution under SAS 
conditions. 

Representative particle differential distributions dN/dp that result from 
the simulation of diffusive acceleration at mildly-rclativistic (internal GRB) 
shocks of speed (3ix = 0.5 are depicted in Figure 1 (see Ellison & Double 
2004, and SBS07) for Fi » 1 simulation results). Here, the subscript x 
denotes components along the shock normal. These distributions are equally 
applicable to electrons or ions, and were generated for ^^scatt ^ 10° , i.e. in 
the SAS regime. Results are displayed for two different upstream fluid frame 
field obliquities, namely 0Bfi = 48.2° and 0Bfi = 59.1°. These define shocks 
with two distinct de Hoffman-Teller (1950; HT) frame dimensionless speeds 
PmT — Pix/ COS ©Bfi ■ The HT frame is the shock rest frame where the flow 
is everywhere parallel to the local magnetic fleld. Note that the distributions 
in the Figure were measured in the normal incidence frame (NIF), the shock 
rest frame in which the upstream fluid flows in along the shock normal; in this 
frame, the magnetic field vectors are generally not parallel to those in the fluid 
or HT frames. The HT flow speed ^iht corresponds to a physical speed when 
it is less than unity, i.e. the upstream field obliquity satisfies cosBba > f^ix ', 
the shock is then called suhluminal. When mathematically /3iht > 1 , no 
Lorentz boost from the local fiuid frame can render the fiow parallel to B and 
de Hoffman-Teller frame does not exist: the shock is said to be superluminal. 

The distributions clearly exhibit an array of indices a , including very flat 
power-laws, that are not monotonic functions of either the fleld obliquity 0Bfi 
nor the key diffusion parameter rj = X/vg. However, it is striking that the 
normalization of the power-laws relative to the low momentum thermal pop- 
ulations is a strongly-declining function of X/vg. This is a consequence of 
a more proliflc convection downstream away from the shock that suppresses 
diffusive injection from thermal energies into the acceleration process. Simu- 
lation runs for X/vg > 10^ inhibit such injection by several to many orders of 
magnitude, and so were not displayed in the plot. Note also that the choice of 
the compression ratio r = 4 is somewhat larger than the Rankinc-Hugoniot 
MHD value for /Six = 0.5 , Mg = 4.04 conditions, and was adopted to afford a 
direct comparison with the semi-analytic convection- diffusion equation results 
of Kirk & Heavens (1989). Details of this comparison are discussed briefly in 
Baring (2009). 
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Fig. 1. Particle distribution functions dN/dp from mildly-relativistic sub-luminal 
shocks {Tixf3ix = 0.577, i.e. f3ix = uix/c = 0.5) of upstream-to-downstream ve- 
locity compression ratio r = uix/u2x ~ 4. Simulation results are depicted for two 
upstream fluid frame magnetic field obliquities, labelled by their corresponding de 
Hoffman- Teller frame upstream flow speeds /3iht = l^ix/ cos 0Bfi . These are in dis- 
tinct groups of three: Qbii = 48.2° (/3iht = 0.75 , multiplied by 100) for the upper 
three histograms, and 0Bfi = 59.1° ( /3iht = 0.975 ) for the lower three histograms. 
Scattering off hydromagnetic turbulence was modeled by randomly deflecting par- 
ticle momenta by an angle within a cone, of half-angle ^scatt , whose axis coincides 
with the particle momentum prior to scattering; three different ratios of the diffusive 
mean free path A to the gyroradius were adopted for each 0Bfi . All results were 
for small angle scattering (SAS), when ^scatt ^ ^/^i and the distributions become 
independent of the choice of ^scatt • A low sonic Mach number Mg was chosen so as 
to maximize the efficiency of injection from thermal energies. 

To summarize the influence of the two key parameters that dictate the range 
of possible spectral indices a , namely the field obliquity 0Bfi , and the ratio 
r] = X/vg , a. parameter survey for diffusive acceleration at (3ix shocks (typical 
of mildly-relativistic systems) is exhibited in Figure 2. Again, the small angle 
scattering limit was employed. To interface more directly with GRB obser- 
vations, it is convenient to represent the particle power-law indices via their 
radiative emission equivalent, namely the high energy spectral index ah (or 
the Band model index f3 = —ah of common usage: see Band et al. 1993), 
for which the observed differential photon spectrum is dn/de^ oc e""'' above 
the MeV-band v-Fy peak. The correspondence between a and ah depends 
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on the emission mechanism and the model assumptions. Here, the popular 
quasi-isotropic synchrotron mechanism (Rees & Meszaros 1992; Tavani 1996; 
Piran 1999; Meszaros 2002) in bursts is adopted for this purpose. Moreover, 
for this illustrative agenda, it is presumed that the emission is in the strongly- 
cooling domain, for which an accelerated electron/pair distribution with index 
a steepens under synchrotron cooling to a power-law of index of a -|- 1 . With 
these specifications, one simply has the relation ct/i = (cr + 2)/2 , so that 
shock acceleration power-laws map over to photon spectra. Accordingly, 
the photon power-law index a^, is plotted as a function of Gefi , with sublu- 
minal shocks constituting those with obliquities ©Bfi < 60° . There is clearly 
a considerable range of possible photon indices spawned by non-thermal par- 
ticles accelerated in mildly relativistic shocks — the same is true for hadronic 
emission components, where, quite often a^'^ u). This is a non- universality 
that can attractively mesh with GRB observations. Note also that the distri- 
butions in Fig. 1 correspond to indices at the obhquities marked in Fig. 2 at 
the top by the labels ^iht = 0.75 and ^mx = 0.975 . 

An obvious feature of this plot is that the dependence of ah on field 
obliquity is non-monotonic. When Xjvg ^ 1 , the value of ayi at first de- 
clines as ©Bfi increases above zero. This leads to very fiat spectra. As ^iht 
approaches and eventually exceeds unity, this trend reverses, and ah then 
rapidly increases with increasing shock obliquity. This dramatic steepening of 
the distribution in near-luminal and superluminal shocks is precipitated by 
inexorable convection of particles away downstream of the shock, character 
that is evinced in Fig. 1. Any amelioration of this rapid decline in the accel- 
eration efficiency requires the reduction of A/r^ to values below around 10. 
Physically, this condition is tantamount to increasing the hydromagnetic tur- 
bulence to high levels that force the particle diffusion to approach isotropy, 
i.e. encroaching upon the Bohm diffusion limit of A/r^ ~ 1 . Then, transport 
of charges across the mean field becomes significant on gyrational timescales, 
and they can be retained near the shock for sufficient times to accelerate and 
generate suitably flat distribution functions. Such low values of A/r^ render 
the flcld direction immaterial, and the shock behaves much like a parallel, 
subluminal shock in terms of its diffusive character. Then, ah is only weakly 
dependent on ©Bfi , another important property illustrated in Fig. 2. 

The third key characteristic is the very flat spectra with ah ~ 1.5 that 
form from extremely flat electron distributions with cr ~ 1 . The origin of 
these (Baring 2009) is in the coherent effect of shock drift acceleration at the 
shock discontinuity. This phenomenon is due to the energy gain of charges 
when they repeatedly encounter u x B drift electric fields (in frames other 
than the HT frame) in gyrations straddling the shock discontinuity; it has 
been widely discussed in the context of non-relativistic astrophysical and he- 
liospheric shocks. When A/r^ ^ 1, the charge trajectories maintain helical 
coherence so that for select gyrophases they can be trapped in the shock layer. 
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Shock Acceleration: Small Angle Scattering 
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Fig. 2. Photon power-law indices ah (for dn/de^ oc £7j°'^ ) corresponding to shock 
acceleration simulation runs in the limit of small angle scattering, again for an up- 
stream flow speed jiix = uix/c = 0.5 , and a compression ratio r = 4 . The indices 
are displayed as functions of the fluid frame field obliquity 0Bfi , with obliquities 
©Bfi > 60° constituting super luminal shocks. The displayed simulation index re- 
sults were obtained for different diffusive mean free paths A parallel to the mean 
field direction, namely A/r^ = 1 (squares), A/r^ = 10 (triangles), A/r^ = 10^ 
(pentagons), and A/r^ = 10'^ (triangles), as labelled. The heavyweight horizontal 
line labelled GRB 080916c indicates the approximate spectral index that is ap- 
propriate for this Fermi burst in the selected time intervals c and e listed in Abdo 
et al. (2009a). All photon indices apply to a cooled synchrotron emission scenario: 
see the text for a discussion. 



gain energy and be reflected upstream to subsequently participate in repeated 
shock encounters. Particles can then efficiently gain energy by the shock drift 
process before they are eventually lost downstream, and the result is a dis- 
tribution function approaching dN/dp oc . Reducing A/r^ and thereby 
introducing extremely modest amounts of turbulence and associated cross- 
field diffusion disrupts this coherence, removes particles from the shock layer, 
and accordingly steepens the spectrum. 
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3 GRB Observations as Probes of Shock Acceleration 



The shock acceleration theory results presented in Figs. 1 and 2 can now be 
interpreted in the light of prompt GRB observations. The focus here is on the 
high energy spectral index of the power-law above the ly-F^, peak. This 
makes direct connection to data from CGRO's EGRET telescope, and now to 
the growing database of Fermi LAT burst detections. To complement these 
probes of diffusive shock acceleration modeling of gamma-ray bursts, it must 
be noted that considerable insights can also be gleaned from consideration of 
the burst spectral shape at and below the u- peak. This was the principal 
focus of the investigation by Baring & Braby (2004), that built upon the early 
work on broad-band spectral fitting of bursts was provided by Tavani (1996). 
The constraints on particle distributions, shock acceleration interpretations 
and viable radiation mechanisms in bright CGRO bursts that were derived in 
Baring & Braby (2004) are reviewed in Baring (2009). 

It should be noted that forging a direct connection between the observed 
afi and the underlying non-thermal particle distribution index a is possible 
when the Thomson optical depth Tt is small, and is straightforward for most 
bursts with good spectroscopy and well-established power-laws above 1 MeV. 
The Fermi-L AT bursts GRB 080916c (Abdo, et al. 2009a) and GRB 090510 
(Abdo et al. 2009b) are suitable examples. Occasionally, the situation may be 
more complicated. For example, some models adopt high Tt (e.g Meszaros 
& Rees 2000; Pe'er & Waxman 2004; Ryde 2005) that lead to strong Comp- 
tonization/thermalization and high compactnesses that spawn rampant pair 
creation. These can prove viable for isolated instances when the observed spec- 
trum is quasi-thermal (e.g. Ryde 2005), and a jet photosphere is invoked. A 
case in point is time interval b in the Ferrni-LAT burst GRB 090902b (Abdo et 
al. 2009c), which offers a highly-peaked and narrow spectral component that 
may (or may not) preclude a shock acceleration interpretation. Yet, this is an 
exceptional selection, and most GRB observations offer fairly clean ah < 3.5 
determinations that permit inferences of an underlying electron or proton non- 
thermal distribution index a . In the majority of extant burst detections, the 
non-thermal spectroscopic interpretation still remains a simpler explanation 
than does the non-isothermal convolution of spectra close to the ideal Planck 
form. 

The most extensive database for high-energy photon spectral indices in 
GRBs is for the CGRO EGRET and Comptel bursts. The EGRET index 
distribution (e.g. Dingus 1995) is constituted by a handful of sources with 
indices scattered in the range 2 < a/^ < 3.7 . Of these, EGRET spark chamber 
detections were concentrated in the index range at < 2.8 , as tabulated in 
Baring (2006). One of the flattest EGRET burst spectra was for GRB 930131, 
with an index ct/j ~ 2 , which corresponds to the indication of X/vg < 10 
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from Fig. 2 (or almost luminal, but not supcrluminal shocks) if a synchrotron 
or inverse Compton cooling model is adopted. The CGRO BATSE at index 
distribution (Preece, et al. 2000) is similarly broad, but with 1.5 < a/i ^ 3.5 
and far greater statistics. The recent Fermi detection (Abdo, et al. 2009a) of 
GRB 080916c in both the GEM and LAT instruments offered a high energy 
index of ah ~ 2.2 (marked in Fig. 2) in its most luminous epochs, and a 
generally steeper spectrum at other times. Accordingly, obscrvationally, shock 
acceleration models must accommodate a radiation spectral index in the range 
2 ^ tt/i ^ 4 in order to be viable. Moreover, they must reasonably account for 
the spectral variability identified in GRB 080916c, i.e. fluctuating ah values. 
This connects to a new development afforded by the refined sensitivity for 
spectroscopy of Fermi: spectral evolutionary studies in temporal sub-intervals 
are now a reality, providing a boon for model diagnostics and a honing of the 
burst paradigm. 

Fig. 2 displays photon indices for electron synchrotron emission in strongly 
cooling scenarios following shock acceleration, where ah = (o"+2)/2 . It is clear 
that the flattest spectral epochs for GRB 080916c are best described by highly- 
oblique, mildly-relativisic shocks, but not quasi-perpendicular ( ©Bfi ^ 70° ) 
superluminal ones. This can also be inferred for the time-integrated spectral 
behavior of the EGRET burst GRB 910503 (sec Baring, 2006). In contrast, 
flat-spectrum EGRET bursts GRB 930131 and GRB 950425 with ah < 2 
cannot correspond to superluminal acceleration regimes if radiative cooling is 
rampant, but demand subluminal shocks with moderate to strong turbulence 
therein, i.e. X/rg ~ 3 - 10 . Other EGRET bursts such as GRB 910601 and 
GRB 990123 (again see Baring 2006 for tabulated indices), and some sub- 
intervals in the Fermi dataset (Abdo, et al. 2009a) for GRB 080916c exhibit 
indices ah > 2.5 . From Fig. 2, these bursts can be consistent with accel- 
eration at mildly superluminal shocks in the 60° - 70° obliquity range, but 
only if the scattering is strong, i.e. X/rg < 3 — 10. No bursts have so far 
evinced extended power-law spectra flatter than ah ~ 1.95 above the the 
z/ - Fj, peak (the possibility of additional high energy components above 100 
MeV in GRB 940217 [Hurley, et al. 1994], GRB 941017 [Gonzalez, et al. 2003] 
and more recently GRB 090902b [Abdo, et al. 2009c], is a separate issue), 
absolving the need for acceleration in shocks with extremely low turbulence, 
i.e. X/rg ~> 10^ regimes. This is fortunate, since such shocks are inherently in- 
efficient accelerators. Moreover, the generation of field turbulence is a natural 
part of dissipation in shocks, so that nearly laminar fields are not expected. 
Field structures that are devoid of fluctuations are never observed in in situ 
magnetometer measurements at hehospheric shocks (e.g. Baring et al. 1997, 
and references therein). Note that moderate fleld obliquities are quite plausi- 
ble in mildly-relativistic internal shocks, contrasting ultra-relativistic external 
GRB shocks, which are necessarily quasi-perpendicular ( 0Bfi ~ 90° ) due to 
the Lorentz transformation of circumburst fields. 
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These inferences are predicated on the presumption of efficient coohng in 
the burst prompt emission zone. This must operate over a significant range of 
electron momenta, since coohng breaks, by an index of Aah = 1/2, are not 
observed above 1 MeV in burst emission. For example, if the spectra span- 
ning the range 1 MeV - 10 GeV in GRB 080916c correspond to a coohng- 
dominated contribution, the acceleration must operate in short, impulsive pe- 
riods, followed by long cooling epochs that permit the electrons to decline in 
momentum by a factor of 10^ or so. Then the cooling epochs should possess 
durations of around 4 orders of magnitude longer than the impulsive acceler- 
ation epochs. This is a significant constraint on cooling scenarios within the 
internal shock model (Baring 2009). The alternative of uncooled synchrotron 
emission yields a photon differential spectral index given by at = {c + l)/2 
above the i^-Fi, peak. Then the injected distribution must have an index a 
higher than that for cooling models by unity, in order to match the burst 
observations. Specifically, the EGRET GRB index range 2 < 0;^ < 3.7 maps 
over to 3 < CT < 6.4 . This is a profound difference in that it pushes the viable 
shock parameter space into the superluminal range, i.e. at higher field obliqui- 
ties, and Bohm-limited diffusion is observationally excluded. In addition, for 
the less popular uncooled hadronic models, the photon and particle indices 
trace each other, i.e. ah — cr , and a similar conclusion that the environment 
is restricted to modestly superluminal oblique shocks is derived. 

It must be emphasized that the spectral index array in Fig. 2 is quite repre- 
sentative of the range and parametric behavior for general mildly-relativistic 
shocks. While the subluminal/super luminal obliquity boundary varies with 

f3ix , the morphology of the index curves, their trends with \/ry, and their 
rapid transition to large indices when crossing to the superluminal domain 
remain qualitatively the same for shock speeds in the range 0.1 < /3ix < 0.9 . 



4 Conclusion 

This paper has explored constraints imposed on the parameter space for 
diffusive acceleration at relativistic shocks by prompt emission observations 
in gamma-ray bursts. The simulation results presented showcase the non- 
universality of the index of non-thermal particles, spawned by a range of 
shock obliquities and the varied character of hydromagnetic turbulence in 
their environs. This non-universality poses no problem for modeling GRB 
high-energy power-law indices ah ■ The observations generally constrain the 
shock parameter space to oblique, subluminal or highly-turbulent superlumi- 
nal shocks not far the Bohm diffusion limit, i.e. X/vg < 10. This presumes 
rapidly cooling synchrotron or inverse Compton emission scenarios. If cool- 
ing arises on timescales exceeding those pertinent for acceleration, then the 
acceleration must occur in mildly superluminal, turbulent shocks; subluminal 
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shocks produce particle distributions too flat to accommodate the observed 
spectral indices. The high photon count detections of the long duration GRB 
080916c and the short burst GRB 090510 in both Fermi's GEM and LAT tele- 
scopes has afforded a new opportunity to constrain shock acceleration param- 
eters in temporal sub-intervals. The prospect of a number of similar platinum 
standard, broad-band GRB detections by Fermi that permit time- dependent 
spectroscopy should hone our understanding of the connection between shock 
acceleration and prompt emission. 
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